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Abstract. The ground state energy of current carriers in graphene considered as a zero- 
gap semiconductor was calculated in the two-band approximation. The condition of the 
electronic (hole) system stability in graphene was obtained. The possibility of the zero-gap 
semiconductor-semimetal transition was discussed. 

It is known that thin graphite films exhibit semimetallic properties [1]; however, a single- 
atomic layer of carbon atoms forming a regular hexagonal lattice (graphene) has such a band 
structure that the energy gap is zero at six K points of the Brillouin zone. Therefore, graphene 
can be considered as a two-dimensional zero-gap semiconductor or a semimetal with zero con- 
duction and valence band overlap [2] . The former approach makes it possible to describe current 
carriers in graphene within the two-band Dirac model^ [7, 8] 

MCK . = ep^, (1) 



where a. = \ are Dirac a-matrices, p = — i/iV is the two-dimensional momentum 




operator (hereafter h= l),u = |7ao = 9.84-10'' cm/s is the quantity similar to the Kane matrix 
element of the interband transition rate, 7 ~ 3 eV is the band parameter, and ao = 1,44 A is 
the interatomic distance in the graphene lattice [9] . In the vicinity of K points of the Brillouin 
zone, the dispersion relation of current carriers is linear, Sp = ±up (+ and — signs correspond 
to electrons and holes, respectively). 



^ratnikovOlpi . ru 

^Dirac equation (1) is equivalent with accuracy of the unitary transformation of the Haniiltonian and the 
wave function of a pair of Weyl equations (see book [3], p. 79). As is known, the Weyl equation describes the 
two-component neutrino in quantum electrodynamics (QED) (see, e.g., book [4]). The use of the Dirac equation 
as a 4 X 4 matrix equation in the two-dimensional system is possible since the 4x4 matrix representation in case 
of two spatial dimensions can be used equivalently with 2x2 matrix representation (see book [5], chap. XIV). 
This fact allows us to extend the formalism of the QED diagram technique to the case of the two-dimensional 
system of Dirac fermions (graphene) . The Weyl equation was first applied to the problem of describing current 
carriers in a zero-gap semiconductor in [6]. 
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For a two-dimensional electron (hole) gas arising during electron doping of a zero-gap semi- 
conductor [2] , the ground state energy per one particle is the sum of three terms 



^gs ^kin ~l~ ^exch ~l~ ^corr^ (2) 



where Ekin = \upF is the average kinetic energy, pp = y ^^"^^ is the Fermi momentum, 
is the two dimensional particle concentration, u is the degeneration multiplicity^. If the 
Fermi level Ep lies above E = 0, the system contains only electrons as current carriers in the 
conduction band with the number of valleys = 2; if Ep < 0, the system contains only holes as 
current carriers with Uh = 2. The Fermi level position can be varied by applying an electric field 
[2]. We can see that both cases in the Dirac model are equivalent. In follows, for definiteness, 
we shall consider the case of electrons. 

The exchange energy is given by the first-order exchange diagram (Fig. 1) 



(Ppde (P'kdijj 
2n2Dj (27r)3 (27r)3 



i?e.ch = -— / -j£^^^Sp{YG{p,e)YG{k,u;)}D^,{p-k,e-u;), (3) 



where the photon propagator D^^, {p — k, e — u) ~ V(p — k) 5^4(5^4 (we neglect the photon pole 
at uj = ±c |p — k|, whose contribution to the integral in frequencies e and u is of the order of 
in comparison with the contribution of Green's function poles), V{q) = , is 
the Coulomb law in the two-dimensional case, and Kg// is the effective permittivity of graphene. 
The free-electron Green's function at A = is [10] 

G(p, e) = - , .""^ + 2m5 {uV - e') N,up, (4) 

where p = (/? = 0, 1, 2) is the convolution with Dirac matrices 7'^ = —i^^o^ for A; = 1, 2 

fl 0\ I 1, X > 0, 

and 7 = , is the 2x2 unit matrix, Np = 9{\p\ — pp)9{e); d{x) = < 

Vo -V I 0, X < 0. 



Expression (3) is transformed to the form (see the Appendix) 

a*Ii /n2D\V2 

^exch = 7=^ , 

V27r V z/ / 



(5) 




k, u 

Fig. 1. First-order exchange diagram. 



^In the general case for two spin components, the degeneration multiphcity is ly — v^^h- It will be shown below 
that the spin-unpolarized state is more energetically favorable than spin-polarized one, for which v — )■ v^j-ijl. 
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where L = f dx f dy f dy ^H'^"'' ^^^^ = = | 

2 

stant [111, a* = — — is the analogue of the fine structure constant. 
The correlation energy is given by the formula [12] 



(G + I) , G = 0, 915965 ... is the Catalan con- 



d'^kdu f dX r -XuV (k) (k, iu) 

2n2D J (27r)3 J T [l - Az/l^ (k) U44 (k, zc^; 




+ Az/V^ (k) Uf^ (k, zcu 



The total polarization operator is written as 



(6) 



(k, iu) = niJ (k, iu) + n^4^ (k, iu) + . 



(1) 



(7) 



which corresponds to the sum of the diagrams 



'T = 1 






The polarization operator Ii^^l{\i,iuo) in the lowest order in the interaction for the two- 
dimensional case is given by [10] 



n£) (k, lu) 



16 



(kp 



,2 |k|"e: 



2.2 
p 



(2vr)^ 



2£. 



(|kr + ^)^-4(k.p-^ 



(8) 



We calculate the correlation energy using a method similar in many respects to the known 
Nozieres-Pines method [13, 14] which is applied to calculate the electron gas and electron-hole 
liquid using asymptotics of the dimensionless polarization operator [15] 



nl?(r,e) 



i-|cose| 



(9) 



where dimensionless variables r = ^/uj"^ + u'^k'^/upp and sin^ = uk/y/to'^ + u^k'^ are introduced. 

To determine the smallness parameter of expansion (7), we estimate the first-order correction 
1144 in the interaction to the polarization operator (8) 



ni'4) (k, lUJ 



d'^pde d'^qd^l 



Sp{rG{p + q,e + n)rG{p,e) 



(27r)3 (27r)3 



(10) 



d^pd^q 2Txa* 
(2vr)4 |q| 



0{pF-\v\)0{pF-\q\) 



a pf 
(87r)2 V 
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Taking into account that the main contribution to 1144'' (k, ion) is made by small transferred 
momenta due to ^(q) [15], (10) should be compared with the asymptotics Ilf^(k,iuj) at small 
|k|, from which we obtain that (10) is small in the parameter ^ 1, which is simultaneously 
the condition of the applicability of the random-phase approximation. 

Substitution of the polarization operator asymptotics at r 3> 1 into (6) yields the contribu- 
tion of large momenta 

3 

^^orr = - ^J^.! [(1 + 9i) ■ ln(i + gi) - gi], (11) 

where gi = ^ ^ 1 at z/ = 2 and a* < 1; expanding (11), we obtain 

^,*2,,l/2 

^corr — „„ . nz-'^"'2D^ \^^) 

8o4v27r 



which coincides with the contribution of the second-order ring diagram (Fig. 2) 

rn(2) ^^'^ f(fpd'^q(PkdedQduj^ r / n 0^ / 1 m 
" J (27r)9 ^ ^ {p-k,e-u)} 



X 



X Sp{rG{q,n)YG{q + k,n + u;)}D^^{k,u)D^,{k,u;)^ (13) 

1 ^^^k^-/,nS(k,..)v(k)^^ 



4n2D J (27r) 

Thus, according to the Nozieres-Pines method [13, 14], when calculating the correlation energy 
at large transferred momenta by formula (6), with an accuracy of the terms of the order of 
gf, the analysis can be restricted to the second-order of the perturbation theory. Apart from 
the second-order ring diagram (Fig. 2), let us also consider the second-order exchange diagram 
(Fig. 3) 

^ ^_i^r iW^dedai^ s^ {r G (p. .) y G (k - p. - .) r . 

4n2D J (27r)y 

■G{p-q-k,e-n-uj)YG{Y>-q,e-n)}D^, (q, fi) Dp, (k, u) . (14) 

Evaluation of integral (14) is very laborious; however, estimations show that, as in the nonrel- 
ativistic case, its contribution is positive and is smaller in magnitude by a factor of ^ than 
(13) [16]. Finally, for the contribution of large transferred momenta to the correlation energy, 
we obtain 

*2 1/2 / ^ \ 

, = -^^^ 1 - nn^l (15) 
864 V 2i/y 2^ ^ ' 

Substitution of the polarization operator asymptotics at r <^ 1 from (9) into (6) yields the 
contribution of small transferred momenta 
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Dan(k,Uj) 



q + k, £ + w 




Fig. 2. Second-order ring diagram. 



7r/2 



3(27r)%2D 



Fig. 3. Second-order exchange diagram. 
The following notations are introduced: 
Pi = p - k, ei=e-uj; P2 = P - q, e2=e-Vt. 



d6{l — cos I 



1 — cos , , ^ 
1-92 — - In ( 1 + 



sin 6 



g2{l - cos I 



+ 



sin 6 



1 — cos 9 

2 In ( 1 + 5(2—7 



52(1 — cos6') 52(1 ~ cos6')^ \^ ' sin^ 
where 52 = ^ 1, and expression (16) can be simplified 



(16) 



^0 

corr 



a 



128V2^ V 4 



2 1 ^^2^. 



(17) 



We can see from (15) and (17) that the order of smallness of the correlation energy and exchange 
energy (5) is a*^ and a*, respectively; thus, we obtain the ground state energy of the electron 
gas as a power series of a*, which we cut off after the terms of the order a*^, 



2V2^ /n2D\V2 a*h fn^oy/^ 



a 



*2 



Stt 25 



1 



64V2^ V 8 27 27u 



u (un2D 



,1/2 



(18) 



It can be seen from (18) that the spin-unpolarized state with the degeneracy multiplicity u = Ue^h 
is more energetically favorable than the spin-polarized state with the degeneracy multiplicity 
= ^e,h/2. The main contribution is given by the kinetic energy; therefore, Egg > 0. The 
condition Egg > in the case of A = means that graphene as a zero-gap semiconductor is 
stable (at Egg < 0, it would be favorable to generate electron-hole pairs). The transition to the 
semimetallic state, i.e., spontaneous electron-hole pair generation, occurs at certain Oq such 
that Egs < at a* > ctg. Equation Egg = upon the substitution of u = 2 yields 



VA^ + B- a = 1.1044, 



(19) 



where A - 3^_^202/27 ^ - 97r-202/9- 
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Let us now consider the effect of some parameters on the band structure of graphene in 
more detail. 

Effective multivalley structure. A structure (superlattice) containing N graphene layers, in 
the absence of transitions between layers, effectively contains v = v ■ N electron (hole) valleys 
with the number of electrons (holes) N e{h) = Ne{h) ■ N, where Ne{h) is the number of electrons 
(holes) in each graphene layer. Let graphene layers be interfaced by a wide-gap semiconductor 
(insulator). 

In the general case the Coulomb law for periodic structures is given by the expression [15] 



27re^ sinh(|q|(i) 



K |q| cosh (|q| d) — cos 



w ' 



(20) 



where < w < 2tt; however, in the case of large transferred momenta such that |q| d ^ 1, the 
second fraction in (20) tends to unity. Formula (11) for -E^^r remains valid with an accuracy 



of the substitution gi 
layers), 



9i 



a u 
12 



corr 



^ 1 and ^2Z) 

1/2 



Ne{h)/ S = N ■ {S is the area of 



36^2^^ ^ I' y 



In 



a u 
~V1 



N 



(21) 



At small transferred momenta, |q| d <^ 1, the Coulomb law (20) is not singular at |q| = 0. In 
this limit, -Ecorr small in comparison with -E^^r '^^^ appearance of In A^; therefore, it can 
be omitted. The ground state energy is given by 



E. 



2V27r 



a 



/i + — In 
36 



a u 



N 



(22) 



Electric field effect. As shown above, application of an electric field can produce a nonzero 
electron (or holes when the electric field direction is changed) density in graphene on substrate 
[2]. In this case, the two-dimensional current-currier concentration is proportional to the gate 
voltage Vg [2] 

= 4^^- (^^) 

where e is the substrate permittivity. 

According to the calculation [17], in graphene containing several layers, the energy gap 
opens at K points of the Brillouin zone under a sufficiently strong electric field. This means 
that the system transits from a semimetallic state (studied in [1]) to the semiconductor state. 
Let us find out whether or not a similar phenomenon exists in single-layer graphene. Assume 
that application of a rather weak electric field results in opening the gap, A 7^ 0, small in 
comparison with the Fermi energy Ep = upp: A ^ Ep = \\/~^^g] in this case, the correction 
to the kinetic energy is = ^ jEp. According to (Al), the correction to the exchange 

energy is given by 

a* Jo A^ 

5i^e.c/^ = -^ — , (24) 



27r Ep' 
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1 1 27r 



where Jo = ( dx ( dy ( d\ , ^ = IQG. Let us express the correction to the correla- 

tion energy in terms of the correction to the polarization operator (k, iu) = U^^^ (k, iuj) + 
51144^ (k, iuj), omitting the terms giving a zero contribution to SEcorr (odd-numbered in u) 

m£) (k,..) = -|^-r^^4^-tan (25) 
fd'kdu VHk)ufhk,zu) ^^(0),,^ . 

dEcorr = / 77T^ [ Cjliy (k, iw) . (26) 



2n2D J (27r)3 i _ (k) nj? (k, lu) 

Changing to dimensionless variables r, 6 and substituting the asymptotics 1144^ (r, 6') at r ^ 1 
and r <^ 1, noticing that the main contribution to 6Ecorr given by small transferred momenta 
due to V^(k) in (26), we obtain 

97ra*V, / 8 \ A2 
5E,orr = ^^In 1 + — ■ — . 27 
256 \ a*^ J Ep 

The correction to the ground state energy contains the additional small factor 

5Eg, = 6Ek,^ + 5E,^,f, + 5E,,rr= [l - — + ^^^J^ In f 1 + ^-ll ^ ~ 0.1202|^; (28) 

|_ TT 256 \ a*^ J \ Ep Ep 

therefore, even in the presence of a gap (small in comparison with Ep), expression (18) is correct; 
moreover, this correction is positive, which suggests that the zero-gap semiconductor phase is 
stable relative to the transition to nonzero-gap semiconductor phase in an external electric field. 

Transition to the semimatallic phase. For graphene, the parameter a* can efficiently vary 
due to image forces, i.e., the variation in the effective permittivity Kg// of graphene depending 

— 1/2 

on its environment (insulator or vacuum); at the substrate thickness / ^ n^pj , it is given 
by Kg// = where e' is the permittivity of a medium above graphene [18]. For a Si02 
substrate, Kg// = 5 and a* ~ 0,44; for a SiC substrate, Kg// = 3 and a* ~ 0,73 [19]. 
Let us estimate the valence and conduction band overlap 5E in the semimetallic state (Fig. 
4). Let initially Ep > and A = 0; and after the transition, the conduction band was 
lowered with respect to the level E = hj 6E/2 and the valence band rose by 6E/2; then the 
number AA^^e of electrons transferred to the valence band is ^PiS; the new Fermi momentum 
is P2 '■ N'^ = ^p^S = Ne + AAi'e, where = ^PpS {pp is the Fermi momentum before 
the transition). We can estimate upi ~ 6E/2. The average kinetic energy of electrons is 
^kin — i'^Pa; for holes, Ej^-^ = ^upi (the number of holes is equal to the number of transferred 
electrons, A^^ = AA^g; therefore, their Fermi momentum is equal to pi). The ground state 
energy in the semimetallic phase is given by 

E'gs = ^kin + Ekin + ^exch + E^orr, (29) 
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K K' K K' 

Fig. 4. Transition to the semimetallic state of graphene depending on a substrate material: the 
conduction band is overlapped with the valence band at all three pairs of neighboring points K 
and K' of the Brillouin zone. 

where Ecorr is calculated by the polarization operator 1144^ = n44-''^ + n44-*''; Eexch is set equal 
to (5) if Nh — Ne so that -Eex-c/i — ^txch — ^exch- Neglcctiug Ecorr, ou the assumption that 
E'gg = 0, we obtain 

5Ec^{h- upF, (30) 

where h = ^^^^"'^^ = ^^^^ = ~ [G + The condition 6 > 1 means that the zero- gap 
semiconductor phase is unstable {Egs < 0); in this case, 6E > 0, i.e., the transition to the 
semimetallic state, the transition from the spin-unpolarized to the spin-polarized occurs. 
This study was supported by the "Dynasty" foundation for noncommercial programs. 

Appendix 

Using the Green's function from [10], we obtain the expression for the exchange energy at 
arbitrary energy gap A 

2n2D J (27r)4 SpSk n^ff \p-k\ 

From (Al), in the nonrelativistic limit, A ^ upF,Bp6]i ^ A^ ^ |p • k|, we obtain 

pnonrel ^ f d^pd'^k 2716^ 



E::T = / ^ iPF -\p\)e ipF - |k|) , (A2) 

n2D J (27r)4 Keff |p - k| 

which coincides with the known expression (see, e.g., [16]). In the ultrarelativistic limit, A ^ 
upF,£p£k ~ IpI ■ |k| ^ A^, we obtain 

K^r- = / ^ f 1 + (p, - iPi) « (P. - iki) . (A3) 



2n2D J (27r)4 ^ |p| ■ |k| / ACe// |p - k 

Dedimensionalizing the integrand in (A3) and integrating over momenta, we obtain the answer 
in the form of (5). Expression (A3) is equivalent to formula (7) in [20]. 
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